Pacotes
library(rgdal)
library(sf)
library(leaflet)
library(plotly)
library(ggplot2)
library(dplyr)
library(stringr)
library(tictoc)
Shape files
municipios_shp <- st_read(dsn = "./Data/rj_municipios", layer = "RJ_MUNICIPIOS_2017_CENSOAGRO") #municípios do estado do RJ
## Reading layer `RJ_MUNICIPIOS_2017_CENSOAGRO' from data source `D:\beatr\Documents\R\CovidEnce\Data\rj_municipios' using driver `ESRI Shapefile'
## Simple feature collection with 92 features and 2 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -44.88932 ymin: -23.36893 xmax: -40.95852 ymax: -20.76321
## geographic CRS: SIRGAS 2000
bairros_shp <- st_read(dsn = "./Data/Limite_de_Bairros-shp", layer ="Limite_de_Bairros")#bairros do município do RJ
## Reading layer `Limite_de_Bairros' from data source `D:\beatr\Documents\R\CovidEnce\Data\Limite_de_Bairros-shp' using driver `ESRI Shapefile'
## Simple feature collection with 163 features and 14 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: 623577.4 ymin: 7446298 xmax: 695396.8 ymax: 7483424
## projected CRS: SAD69 / UTM zone 23S
limite_muni <- read_sf('./Data/Limite_do_Municipio_do_Rio_de_Janeiro-shp/Limite_do_Municipio_do_Rio_de_Janeiro-shp.shp') #limite do município do RJ
limite_uf <- read_sf('./Data/RJ/33UFE250GC_SIR.shp') #limite do estado do RJ
Mapa dos municípios do Rio de Janeiro
tic()
ggplot() +
geom_sf(data = municipios_shp, colour = "dark blue", fill = NA) +
geom_sf(data = limite_uf, colour = "dark blue", fill = NA, lwd = 1) +
theme_void() +
ggtitle("Mapa do estado do Rio de Janeiro")
toc()
## 1.39 sec elapsed
Mapa dos bairros do Rio de Janeiro
tic()
ggplot() +
geom_sf(data = bairros_shp, colour = "dark blue", fill = NA) +
geom_sf(data = limite_muni, colour = "dark blue", fill = NA, lwd = 1) +
theme_void() +
ggtitle("Mapa do município do Rio de Janeiro")
toc()
## 0.97 sec elapsed
Gerandos dados aleatórios
dados_bairros <- data.frame(NOME = bairros_shp$NOME, Valor = round(rgamma(dim(bairros_shp)[1], shape = 2, rate = 0.5),2))
dados_municipios <- data.frame(NM_MUN = municipios_shp$NM_MUN, Valor = round(rgamma(dim(municipios_shp)[1], shape = 7.5, rate = 1),2))
Plotando os mapas com os dados gerados
theme_map <- function(){theme(
panel.ontop = TRUE,
panel.grid = element_blank(),
line = element_blank(),
rect = element_blank(),
text = element_blank(),
plot.background = element_blank())}
tic()
ggplotly(municipios_shp %>%
left_join(dados_municipios, by = "NM_MUN") %>%
ggplot() +
geom_sf(aes(text=paste('</br>Município: ', str_to_title(NM_MUN)), fill = Valor), color="dark blue", show.legend = FALSE) +
geom_sf(data = limite_uf, colour = "dark blue", fill = NA, lwd = 1) +
theme_map() +
ggtitle("Mapa dos municípios do Rio de Janeiro"))